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Abstract. 

A comprehensive input-output theory is developed for Fermionic input fields. Quantum 
stochastic differential equations are developed in both the Ito and Stratonovich forms. 
The major technical issue is the development of a formalism which takes account of 
anticommutation relations between the Fermionic driving field and those system operators 
which can change the number of Fermions within the system. 



1. Introduction 

The introduction of Input-Output formulations in the 1980s [l][3]|4] was a response to the 
necessity for a theory of quantum damping which could deal with travelling wave situations. 
A formulation of a photodetector theory was suggested using input-output methods EH5J, 
in which the detection of a photon is envisaged as the conversion of an input light field into 
an output electron field. However, to do this requires a theory of inputs and outputs with 
Fermion fields, and an elementary theory was developed. Although basically satisfactory, this 
formulation was not complete, and in particular, could not be used to derive a master equation. 
The problem that arises is quite simple: the equations of motion for system generators are 
different, depending on whether a system operator is viewed as commuting or anticommuting 
with the fermionic heat bath operators. 

In this paper this technical problem is overcome. It is shown how we may define 
"restricted" system operators so as to commute with all bath operators, and that these internal 
system operators all obey a quantum Langevin equation of the same form. 

However, the equations of motion in the original operators are, for the two level atom, 
linear and thus exactly soluble. Thus, we have a description of a two level system interacting 
with a Fermi Bath which is essentially the same as that of a harmonic oscillator interacting 
with a Bosonic heat bath. The two level atom behaves like a Fermion coming to equilibrium 
with all the other Fermions. 

It is possible to develop Fermionic quantum stochastic integration, and the corresponding 
Ito and Stratonovich formulations of quantum stochastic differential equations, and finally, 
from these, to derive the Master equation in the expected form. 

To the best of the author's knowledge there have only been two other papers dealing 
with input-output theory of Fermions, that of Sun and Milburn 1 8 1 and that of Search et al 
171 . Neither of these attempts a comprehensive input-output formalism; the first concentrates 
essentially on counting statistics, while the focus of the second paper is on the theory of 
Fermions inside and outside of a linear cavity, and does not address the quantum stochastic 
issues which are the principal topic of this paper. 
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2. Beams of Fermions 

2.1. Input and Output Fields 

We want to consider here non-relativistic Fermi fields, in which we shall for simplicity and 
clarity make no mention of spin (though this is always present in Fermion, it plays no essential 
role in their non-relativistic description). For simplicity, we shall also consider propagation in 
only one dimension. A Fermi field can then be written: 

(1) 



/ i r 

D(t,x) = \ — dke^-^dik) 

V 27T J_ QO 



where d(k) is the destruction operator, with anticommutation relations 

[d(k),rf(k')] + = S(k-k') 

[d(k),d(k')} + =[<fi(k),<0(tf)]+ = O (2) 
leading to the equal time anticommutation relation, 

[D(t,x),D^(t,x')] + = S{x-x'). (3) 
The dispersion relation 

uo = hk 2 /2m (4) 

follows from non-relativistic mechanics, and ensures D(t,x) obeys the time dependent 
Schrtidinger equation 

m 9D(t t x) = h 2 d 2 D(t,x) 
dt 2m dx 2 

In the situation we wish to consider an input field radiating into a system, which itself radiates 
an output field in the opposite direction. In that case, it is more appropriate to consider the 
field as being defined on the half-line < x < oo, for which an appropriate expansion is 

D{t,x) = <J^J^ dkcos{kx)e iult d(k) (6) 
which can also be written 

/ 1 /-oo 

D(t, x) = \ — dk{e^ kx -^d in (k) + e~^ kx+ ^d out (k)} (7) 

V ^7T J Q 

where 

dta(fc) = dout(fc) = d ( k )- (8) 
We can also define 

Din(t,x) =J— dke^ kx -^d in (k) 

V 2ir J Q 

D out (t,x) = J— dkc- i( - kx+ ^d out (k) (9) 

V 2ir J 

and the boundary condition 

D in (t,0)-D out (t,0) = 0. (10) 

follows from J1J. 
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2.2. Dispersion 

Matter waves are dispersive; in fact the group and phase velocities differ by a factor of 2 for all 
frequencies. This means that the simple propagation of statistics as in light beams is not valid 
here — there is no solution analogous to the solution of the one dimensional wave equation 
which can be written in the form 

A(t,x)=A in (t + x/v). (11) 

We will therefore limit our considerations to very narrow bandwidth situations. In order to 
judge what sort of bandwidth can be considered "narrow", let us consider the correlation 
functions. 



2.3. Correlation Functions of Propagating Fermion Beams 

We want to consider correlation functions 

G^(x 1 ,t 1 ;x 2 ,t 2 ) = (D\t 1 ,x 1 )D(t 2 ,x 2 )) (12) 

and 

G {2 \x 1 ,t 1] x 2 ,t 2 ;x 3 ,t 3 ;x 4 ,U) = (D\x 1 ,t 1 )D^(x 2 ,t 2 )D(x 3 ,t 3 )D(x 4 ,t i )) (13) 

These are analogous to the similarly defined correlation functions for optical fields. 
Let us now consider a stationary narrow bandwidth field, such that 

(a f (k)a(k')} = N(k)6(k - k') (14) 



where 

Now defining 
then 



N{k) = unless k w fc . (15) 

N(u)dw = N(k)dk (16) 

(DUt,x)D(t',x')) = — [ duNl^e^-^-^-^W (17) 
2n J 

If u>o, ko are the central frequency and wavenumber of the range, and the central velocity of 
propagation is 

v = uja/ko (18) 
then, provided that the range of frequencies, Su>, satisfies 

(x - x') 5u/v < 1 (19) 

we can first write 

g(T)=e- i »° T (Dt(T,x)D(0 > x)) (20) 
and then we can derive 



(21) 



G {1) {x,t;x' ,t') = e -o[t-t'-(*-*')H f x + ^^JQ g(t-t'- ^-^) 

This equation gives the correction due to dispersion of the "propagation approximation" to 
evaluating the correlation function of a fermion beam at two separated points. The condition 
of validity for such an approximation is dl9> , which can be interpreted to mean that the size 
of individual wavepackets, which is order of magnitude v /8cu, must be very much larger than 
the distance \x — x'\ between the points considered. 

The correction in (12 It is relevant to the measurement of time correlation functions by 
delayed coincidence measurements, when the delay is induced by allowing one beam to 
propagate further than the other. 
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2.4. Thermal correlation functions for Fermionic beams 

In the optical case a thermal light beam is considered to be Gaussian, and the factorizable 
property of Gaussian moments leads to a relationship between 1st and 2nd order correlation 
functions which is the characteristic of the "bunched" nature of thermal light. It is difficult 
to define what might be considered to be a "Gaussian" Fermion state, but a thermal Fermion 
state can be defined. 

In this case of a thermal state we have 

(d\k)d{k')) = 5{k-k')N{k) (22) 

(d}{k)d ] (k')d{k")d{k"')) = [5(k - k")5(k' - k'") - 5{k - k'")5(k' - k")] N(k)N{k'). 

(23) 



Here, if 



then 



N(uj) duj = N(k)dk, (24) 



N(u) = l/ [cxp{huj /kT) + 1] . (25) 

The antisymmetric requirement (I23> is in fact the natural analogy of the Gaussian factorization 
property of photon beams (there would be a + sign on the RHS of (I23> in the case of photons 
instead of a — sign.) It leads to the relationship 

G (2) (xi, h; x 2 , t 2 \ x 3 , t 3 ; x 4 , U) = G (1) (zi, h; x 3 , t 3 )G (1) (x 2 ,t 2 ; x 4 , U) 

- G« (zi, h; Xi , t 4 )G (1 ' (x 2 ,t 2 ;x 3 ,t 3 ). (26) 

The corresponding formula for a Gaussian Boson beam differs only by having a positive sign 
rather than a negative sign on the right hand side. 

In the case of stationary statistics, evaluated with time difference t at x = 0, we have for 

3 (2) (r) = G (2) (0,i;0,t + r; 0, i + r; 0, 0) (27) 

g ^( T ) = G (1) (0,i;0,t + r) (28) 

that 

9 {2 Hr)=9 {1) (r)9 {1 \-T)-[g^(0)r (29) 

= l5 (1) (r)| 2 -| 5 ( 1 )(0)| 2 (30) 

Clearly .g' 2 '(0) = 0, corresponding to perfect antibunching, as expected from a Fermion 
beam. 



3. Interaction of a System with a Fermionic heat bath 

In order to fix our ideas, let us consider the ionization of an atom under the influence of an 
impingent electron beam. The Hamiltonian is 

H = H sys + Hi n t + Hb ■ (31) 

Here, H sys is the free atom Hamiltonian, whose precise form will be left open. The bath 
Hamiltonian, corresponding to a field on a half line, < x < oo, can be written 

/>oo 

H B = dkl%oj{k)d\k)d{k) (32) 
Jo 

where in this case, for an electron of mass m, 

uj(k) = hk 2 /2m. (33) 
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The Fermion operators d(k) are as in the previous section. Finally, the interaction is conceived 
as representing the absorption or emission of an electron, and is written 



Jo 

Here c, c) are system operators. The action of & on a system state increases the number of 
constituent electrons by 1 . 

3.1. Fermionic and Bosonic System Operators 

The commutation relations for the system operators depend on the systems being studied. 
In the usual case of a Boson bath [1] the system and bath operators commute at equal times. 
However, when dealing with a Fermionic bath the situation is different, depending on whether 
a system operator can be considered as changing the number of Fermions which make up the 
system or not. The separation into "system" and "bath" tends to obscure the fact that the 
system does have an internal structure, and that (for example) an ion and a neutral atom must 
have different numbers of constituent electrons. An operator such as c, which can be regarded 
as removing an electron from the neutral atom, must anticommute with all bath operators, 
since it must be composed of an odd number of creation and destruction operators. On the 
other hand, there are operators (such as H sys ) which do not change the number of constituent 
electrons, and hence commute with the bath operators at equal times. We thus conclude that 
these are two kinds of system operators. 

a) Bosonic — these commute with all bath operators, d(k), d^(k) (at equal times) 

b) Fermionic — these anticommute with all bath operators d(k), d^(k) (at equal times) 

We will use the notation a, b, c, etc. for the system operators to emphasize that such operators 
may anticommute with the bath operators, and hence are not necessarily independent of them. 
We will shortly introduce "restricted" system operators, which are independent of and hence 
commute with the bath operators. 

3.2. Derivation of Quantum Langevin Equations 

Because all Fermion fields of interest are massive, and therefore the wave propagation is 
dispersive it is not quite as easy to derive quantum Langevin equations as in the optical case. 
Added to this is the complication that some of the system operators are Fermionic, and others 
Bosonic. 

The operator c which occurs in the interaction Hamiltonian must be Fermionic, since it 
changes the number of constituent electrons in the system by 1 . Using this, we derive the 
equation of motion for d(k). 




(34) 




iuj{k)d{k) — kc 



(35) 




■i 



d(k,t) 



(36) 



to 



We now define 




(37) 



(38) 
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The quantity d(t) is a genuine Heisenberg operator for the time t, whereas d- m (t) depends 
only on the initial values d(k, to), of the destruction operators. We now write the equation of 
motion for an arbitrary system operator a, 



a= - — [5, H sys ] 



- [a, H sys ] 



dkn(k) [a,d\k)c- c\k)d{k)] 



(39) 



dkK(k){Td^(k)[&,£}± - [a,&]±d(k)} (40) 
where the top signs apply if a is Fermionic, and the bottom signs apply if a is Bosonic. 



3.2.1. The white noise approximation To obtain Langevin equations we must make 
approximations. There are two principal approximations. 

a) The interaction is weak, and the free motion of c(t) is proportional to e~ 1Wot . 

b) The frequency u>o is rather large. Since the equation (I40i is homogeneous in a, this 
means that the main contribution from the integrals will occur where co(k) w u>o- 

We can thus write an approximate expression for the k integrals in d40l by evaluating the 
c-number n(k) at fco; i.e., 

/>oo 

(41) 



(42) 



dk n(k)d{k) « K(k ) / dk d(k) sa V2tt re(fc ) d(t). 
Jo 

provided n(k) is a smooth function of k around ui(k) = ujq. From d36l 

/ 1 r°o 

d{t) = d in (t) + \ — dk' n^e-WW*-* >ctf)dtf 

V 27T ,/n 



d ia (t) 



1 dt' [ ^^M^fc^c-^*-''^^) 
du) 



(43) 

and provided K,(k(u>)) and dk{io)/dio are smooth around u> = u>o, we can again approximate: 



d(t) = d in + V2tt / S(t - t')c(t') 



dk(u>) 
dw 



that is 



dW=*n(*) + y|«(*M)^^2(t) 
Before finally substituting to derive the quantum Langevin equations, notice that 



1 

27 
1 

27 



and if this is being used mostly at frequency ujq ^> 0, then we can approximate to get 



din(t),dl(t') 



dk(u>a) 
dk 



5{t-t') 



(44) 
(45) 

(46) 
(47) 

(48) 
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We want a noise input with anticommutator normalized to S(t — t'); we therefore define 

i 

dk (o>o) 



fin(t) 

so that ( 1411 becomes 



din(*), (49) 



dkn{k)d{k) « |c+V7/in(t), (50) 



where 



7 = 27r K (fc(^)) 2 (^^^J . (51) 

Using this approximation, the interaction Hamiltonian can also be approximated by 

ffmt « i* { (|c + V7/iL(*)) 2 - 2 f + V7/inW) } , (52) 
and this is a form which will be useful in the remainder of this paper. 



/ dk((jja) \ 



3.2.2. Fermionic quantum Langevin equations Now substitute into J40i to get 

«= -^[a, J ff 8ys ]-[a,ct] ± {|c + V7/in(t)}T{|£ t + V7/iL(*)}[«>5]± (53) 

Equations J53i are the Fermionic quantum Langevin equations for the full system operators. 
At this stage all we know is that the fi n (t) are determined by the d(k, to), whose statistics are 
determined from the initial state of the Fermionic heat bath, and that 

[fin(t),fL(t')} + = Kt-t'). (54) 
The validity of (1531 and 1541 is restricted to situations in which the interaction is rather weak, 
and in which the time dependence of c(t) in the case of no interaction is e _IWot , for some 
rather large ui n . 

3.2.3. "Out" operators As in the case of Bosonic quantum white noise, we can define "out" 
operators, by considering solutions of (1351 in terms of a final condition at a time ti > t. 
Thus, ( 1361 becomes 

rti 

d(k, t) = e-^lH)^, h) + J K (fc)e- iw «(*-* h(t')dt'. (55) 



and we define 



/ 1 poo 

d out (t) = J— / dke-WW-^dfah). 

V 27T ,/n 



(56) 



and from d56l and d43l 

d ou t(t) - d in (t) 



-i pti poo 

- dt' dkKiky^^-^cit^dt' (57) 
2tt J to J a 

and using the same methods as in ( 1441 to (1451 

dout(i) - d in (t) = V2^n(k(uj )) ^f^-c(t). (58) 

auiQ 

or, in terms of f in (t), / ou t(*), 

/out(i)-/in(t) = V7c(t). (59) 
We can derive "time reversed" quantum Langevin equations in terms of these "out" noises. 
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3.2.4. Commutation Relations between System and Input Operators The operators 
d(t),d'(t), defined by ( I37i . commute (anticommute) with all Bosonic (Fermionic) system 
operators at the same time, since these describe independent degrees of freedom. Thus, using 
( 145 \ and i49i . we can say that if a(t) is an arbitrary system operator 

[a{t),d(t)] ± = [a(t),f ia (t)] ± = -l^f[a(t),c(t)} ± (60) 



If we rewrite d59l as 



/ln(*) + \Vlc(i) = f out (t) - l -^c{t) (61) 



it is easy to rewrite the quantum Langevin equation (1531 in the "out" form 

^[a,H sys ] - [a,5 f ] ± {~2 £+ \/7/out(*)} 



+VlfL(t)}[a,£}±. (62) 



3.3. Restricted System Operators 



Because the bath operators do not commute with those system operators a which are 
Fermionic, we have different forms for the quantum Langevin equation depending on whether 
or not the operator under consideration is Fermionic. 

We shall introduce a different set of system operators, called restricted system operators, 
which do not have this problem, that is, the equations of motion take the same form for all 
restricted system operators. 

To do this we introduce the operator I, in the bath space, which anticommutes with all 
bath operators. 

[I,d(k)]+ = [I ) <fl(k)]+ = (63) 

This operator is easy to construct explicitly. If \A) is any bath state with a definite number ua 
of bath Fermions, then we define 

7|A) = (-1)^L4>. (64) 

Clearly I is Hermitian, I 2 = 1, and I commutes with all system operators. 
We now define restricted system operators x by 

Ix if x is Fermionic 

x if £ is Bosonic 

Independently of whether x is Fermionic of Bosonic, x commutes with all bath operators. 
These restricted system operators, x, are to be contrasted with the full system operators i. 
The x are the more physical operators, since they are the ones that turn up in the Hamiltonian. 
The restricted system operators are necessary when we wish to consider operations on the 
reduced density operator p sys = Tib {p}- 

In order to rewrite the Hamiltonian in terms of the restricted operators, it is necessary 
only to rewrite Hi nt as 

H Int = ih dkK(k){rf(k)Ic-Jld(k)} (66) 
Jo 
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3.3.7. Commutation Relations for I(t) Because / does not commute with Hi Dt , it is time 
dependent, and we re-write is as I(t). Because I(t) anticommutes with d(k,t), d^{k,t), we 
deduce from ( 145 \ and d47t that 

[/(*), / in (i)] + = - fyc(t) = -yftim (67) 
[I(t), /+ (t)] + = - V7c f (*) = -V7^ f W (68) 
The commutators of I(t) with either restricted system operators x(t) or the full system 
operators x(t) are zero. The fact that I(t) does not commute with f- m {t) or f? n (t), and is 
time dependent, makes the quantum Langevin equations for the restricted operators, 1721 . in 
general quite different from those for the full operators, <53l . 

However, because I(t) 2 = 1, and because I(t) commutes with system operators of either 
kind, all equal times algebraic relations between different system operators are the same for 
both full and restricted system operators. 

3.3.2. Langevin equation for I(t) Using the commutation relations J67I68L the equation of 
motion for I(t) can be deduced by much the same reasoning as above, to be 

t = - 2^f[fl(t)c + C t/ in (t) + ^ylJc) i (69) 

= - 2V7 {fL(t) J £ + &Ifin(]k) + v/t"^} ■ (70) 

Finally, notice that I(t) 2 = 1 implies that II + II = 0. This can be explicitly shown from 

by using i67l68l . 

The initial condition, that is the operator I(to), is determined by the formulae (163164b 

using d(k,to), cft(k,to), and by definition f- m (t), fi n (t), are linear functions of d(k,to), 

d'(k, to). Thus, we can say that 

[I(t o )J ia (t)} + = [I(t o ),fl(t)] + =0. (71) 

3.3.3. Langevin equations for the restricted system operators Similar reasoning can be 
used to deduce that for all restricted system operators a, the quantum Langevin equations 
take the form 

a = -^[a,H sys ] - [o, c t] {Zc+^J/^)} + |2 c t + ^y/t,(i)/} [a,c]. (72) 

This equation can also be deduced by substituting a = a, or la as the case may be, into J53i . 
and using J70i . together with J67t . 

The major advantage of the use of restricted system operators is that they are truly 
independent of the bath operators evaluated at the same time. The whole burden of 
antisymmetrization required between bath and system is borne by the operator I(t), which is 
not purely a function of fi n (t), // n (i), but is a dynamical variable whose equation of motion 
must be considered alongside that of the system operators a — namely the restricted operator 
quantum Langevin equation (I72> . 

3.3.4. Fermion conservation superselection rule Fermions can only be created and 
destroyed in pairs. Hence if Nb is the number of Fermions in the bath, and iV sys is the 
number contained in the system, the quantity 

K=(-1) Nb+n ^ (73) 
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is a conserved quantity. Furthermore, the quantity 

J(t) EE (-1)^ (74) 

is a Bosonic system variable, which therefore possesses the properties 

I(t)J(t) = J{t)I(t) = K, (75) 
[g(t),J(t)}± =0 (76) 



In the latter equation g(t) is an arbitrary restricted system operator, and the choice of 
anticommutator or commutator depends on whether the corresponding full system operator 
operator g(t) is Fermionic (+) or Bosonic (— ). 

This relationship means that correlation functions involving Fermionic full system 
operators can be evaluated using only the restricted system operators and the conserved 
operator K. Thus for example, noting that c(t) is Fermionic, 



(c\t)c(t')) = (I(t)c\t)I(t')c{t')) (77) 

= < C t(t)7(i)7(f)c(f)) (78) 

= (c^t)J(t)K 2 J(t')c(t')) (79) 

= (&{tW)). (80) 
where we have defined the notation 

c=Jc, (81) 



which amounts in practice to a restricted operator expression of the full system operator c. 
Only if there are an odd number of Fermionic system operators in the correlation function 
does the nature of K come into play, and then we find the relation between the c correlation 
functions and those involving c involves a plus or minus sign, depending on whether the total 
number of Fermions in the system is odd or even. Such correlation functions are therefore not 
likely to be of much physical relevance, and in fact involve interfering states with different 
total numbers of Fermions; a violation of the superselection rule. 

3.3.5. Properties of solutions The solutions of ( 1691721 at time t depend on the initial 
conditions at time to, and /i n (f ), /i n (f ) for t < t' < t. In the same way as we derived 
( 17 U . we can show that for all a(t), 

[a(to)J ia (t)} = [a(to),fUt)}=0. (82) 
From the Langevin equations <69l72t we can show that 

[a(i'), /in (*)] =[a(t') 1 fl(t)}=0 forio <t'<t; (83) 

[I(t'),f- m (t)]+ = [I(t'),fl(t)} + = for t <t'<t. (84) 

The proof is quite simple; assume the results are true for some value s of t'. Using the 
Langevin equations evaluate the derivative with respect to s of the left hand sides of d83l84i . 
and then use the relations ( I83l84t for t' = s to show that the results vanish. Hence that the 
only solution of the resultant differential equations corresponds to the truth of ( I83I84I . 
In a similar way, one can show that 

[a(t'),f ont (t)] =KO>/out(*)]=0 fort<f'<t i; (85) 

[7(f), jout(i)]+ = [7(f), / t ut (*)]+ = for t < f < h. (86) 



Input and output in damped quantum systems III: Fermionic. 



11 



Noting now ( I59> . we can derive 

[a(t'),f in (t)] = -u(f -t)^y[a(t'),c(t)], (87) 

Ht'),fL(t)} = -u(t'-t)Sy[a(t'),c*(t)}, (88) 

W), /*(*)] + = - u(f - t)^[I{t'), 5(f)], (89) 

[I(t'),fUt)}+= -u{t' -t)^[I{t'),~c\t)]. (90) 
From these it also follows that for any Fermionic full system operator g 

W), /in(*)]+ = - - t)^[g(f),c(t)}+, (91) 

[ff(*') J /i r n(*)]+ = t)Vim'), 5 f (*)]+• (92) 



These results are derived by writing g(t) = g(t)I(t), and expanding the commutator using 
J87H90I I. The corresponding results for a Bosonic full system operator, which is identical with 
its corresponding restricted form, have the same form as d87l88i . 



4. Fermionic Quantum White Noise and Quantum Stochastic Differential Equations 

The operators /; n (i), /inM> have the idealized anticommutation relations j54t . which leads 
naturally to a formulation of Fermionic quantum white noise. We define a Fermionic Quantum 
Wiener Process by 



F(t,t ) = / hn(t')dt', (93) 

and we assume the averages, 

(F\t,t )F(t,to)=N(t-t ) (94) 
{F(t,t )F^(t,t ) = (l-JV)(f-to) (95) 

and the anticommutator, from d54l 

([F(t,to),Fi(t,t )}+) = (t-t ) (96) 

We of course also assume the independence of the F operators defined on non-overlapping 
time intervals; 

(Fi(t,t )F(s,s )) = (F(s,s )F\t,t )) 

= [F(t,to),F(8,8o)]+ = (97) 
if (s, s ) and (f , i ) are disjoint. In frequency space, this can be obtained by writing 

rr~ r°° 

/*(*) = \hr foHc-^du; (98) 



2tt _ 
with 

[/o(w),/ V)]+ = <*(<<> (99) 
[/iM/om =JV%-w') (100) 
(/oM/ V)> =(1-^)%-^') (101) 
The /o(w) are therefore like idealized Fermion destruction operators, defined on a frequency 
range (—00,00). In practice, the do(k) are the true destruction operators, and the 
correspondence between fo(u>) and do(k) is made via the relationship (1491 . and is only valid 
for a narrow bandwidth around the frequency lu^. The components of fi a (t) and di n (t) outside 
this narrow bandwidth have little effect on the solutions of the quantum Langevin equations. 
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Corresponding to the formula for Bosonic white noise, the density operator which gives 
the averages J94I95I has the form 

fiF^t,t )F{t,t y 



Pw(t,t ) = (1 + e M )exp 



t-t 



(102) 



in which 



N = — — (103) 
+ 1 

The formulation of F(t, to) as above enables us to develop a formal theory of quantum 
stochastic integration and quantum stochastic differential equations which is quite simple and 
easy to use, and whose use gives essentially the same results as any exact formulation. We 
will also be able to show that the quantum stochastic differential equations so developed are 
exactly equivalent to the master equation for systems interacting with a Fermionic heat bath. 
To do this, we partition the time interval (t ,tf), inside which we are interested in treating 
the motion, into subintervals at bounded by times t , 1 1 , t 2 , ■ ■ ■ , t n = tf corresponding to the 
increments 



At, =t i+1 -ti (104) 
AFi = F(t i+1 ,ti)= / + f in (t')dt'. (105) 



Depending on whether g(t) is a Bosonic or a Fermionic operator, AFi will commute or 
anticommute with g(ti). 

The corresponding joint density operator for the Fermi noises partitioned this way is then 
the direct product 

PF = pw{t n ,t n -i) ® (g> ...(g) pw{t 2 ,ti) ® p W (ti,t ). (106) 

This form means that we can consider a trace operation over the Fermion bath to be taken over 
each time interval in a discretization; thus if Q is some operator which acts on the fermion 
bath, we can write 

Tr B {Qp F } = Tr (t „ i t„_ l) {Tr (tn _ liU _ a) {. . . Tr (t2itl) {Tr (tl , to} {Qp F }}}} . (107) 
This form is particularly useful in deriving correlation function identities. 



4.1. Quantum Stochastic Integration 

As in the case of classical stochastic integration with respect to white noise, there are two 
natural definitions of integration, the Ito and Stratonovich methods. The definition of these is 
relatively straightforward. 



4.7.7. Fermionic quantum Ito integral If g(t) is a full system operator (not a restricted 
system operator), the Fermionic quantum Ito integral is defined by 

(I) / g(t')dF(t')= lim VsfeJAf} (108) 
J t V 

where t < tx < t 2 < ■ ■ ■ < t n = t, and the limit is a mean-square limit. A similar definition 
is used for J t g(t') dF^(t'). The advantage of the Ito definition J 1 08l > of the integral is that 
the increment, AFi is seen in the explicit definition on the right hand side, to be in the future 
of U. 
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In particular, as is the case for classical and Bosonic stochastic integration, 

(I) jf g(t')dF(t')^ =0, (109) 

CDjT g(t')dF\t')j =o. (HO) 

Further, depending on whether gr(t) is Bosonic or Fermionic 

(I) f [g(t'),dF(t')] ± =0, (111) 

Jt 

(I) / [ff(i'),rfi ;,t (i')] ± =0. (112) 



4.7.2. Fermionic quantum Stratonovich integral The Stratonovich integral can be defined 
in the same way as it is for Bosonic noise 

(S) f 9 {t>)dF(t') = Urn Yl m \ m+l) Afl. (H3) 

"'to j 

As in the case of classical and Bosonic noise, we cannot make any connection between these 
two forms of integral without knowing what kind of stochastic differential equation is obeyed 
by the system operators. 

Simple relations of the kind ( 1109m 121 do not hold. For the first two we need to establish 
the relationship between the two kinds of integral first, and this will be done in Sect l4.3l 

4.2. Ito quantum stochastic differential equation 

We will define the Ito quantum stochastic differential equation obeyed by a restricted system 
operator a as 

(I) da = --[a,H By8 \ dt + ^(1 - N) {2c f ac - ac< c - c< ca} dt 
+ {2mc t - acc) - cc) a) dt 

- V7 [a, ct] I dF(t) + dF^ (t) I [a, c] . (114) 

This can be written in the alternative form 

(I) da = - ~ [a, H sys ] dt + j -c* [a, c] - ^ [a, c^cj 

+ ^v{[c, [o,ct]] + -[ c t,[o,c]] + } dt 

- [a, c^] I dF(t) + y/jdF^(t) I [a, c]. (115) 

These definitions are made with foresight; using them we will show how to connect the 
Ito and Stratonovich definitions of the quantum stochastic differential equation, ultimately 
showing that the Fermionic Langevin equations as derived above are then to be interpreted as 
Stratonovich quantum stochastic differential equations. 
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4.3. Connection between Ito and Stratonovich integrals 

We write the term g(ti+i) in the definition of the Stratonovich integral as 

g(t i+1 ) =g(k) + Ag{U), (116) 

where Ag(tj) is to be calculated using the Ito quantum stochastic differential equation il 14> . 

Since Ag(ti) is to be used in J 1 1 31 . only the terms involving the noise operators will be 
significant, since in the limit terms involving Ati as well as a noise operator vanish. 

The conversion formula will differ, depending on whether g is Fermionic or Bosonic; we 
will consider first the Fermionic case. The quantum stochastic differential equation il 141 ) is 
written for the restricted system operators, but we can write 

g{t i+1 ) = I{t i+1 )g(t i+1 ) (117) 

= K J(t i+1 )g(t i+1 ) (118) 

= Kg(t i+1 ). (119) 

Since g is a restricted system operator, it obeys the quantum stochastic differential equation 
il 14i . which means that we can write for the stochastic part of Ag(U) 

A.9fe )Locha S tic =K{-^[g,c^IAF t + ^AF^I[g,c}} (120) 

= K{-^flJg,J}IAFi + ^/AF}l[Jg,c]} (121) 

= if {-V7[ff,ct] + J/AJi + ^AF^IJ[g,c] + } (122) 

= K{-Sy\g,c*] + KAFi + AF}K[g,c} + ] (123) 

= K 2 {-^[g^^+AFi - ^jAF}[g,c] + ) (124) 

= - V7[ff,c+] + AJi- ^AF}[g,c] + (125) 
From ( 1941951 we can write in stochastic integrals 

AF}AF t = NAU, (126) 

AFiAFf = (1 - N)At h (127) 
and the anticommutation relations also give 

AFiAFi = AF?AF} = 0. (128) 
Carrying out similar reasoning for other integrals, we find that we can write 

Fermionic integrand g: 

(S) f g(t')dF{t') = (D ( g( t ')dF{t')-^- [\g(t'),c(t')] + dt', (129) 



JtQ Jto Jto 

(S) f dF(t')g{t>) =0) f dF{t')~g{t')- ^ {1 ~ N) f W),c(t')} + dt> , (130) 

Jto Jto Jto 

(S) f ~g{t')dF\t') = (Y) [ ~g{t')dF\t')- ^ {1 ~ N) [ [g(t'),J(t')] + dt', (131) 

Jto Jto z Jt 

(S) f dF\t')~g{t') = <\) f dF\t')g{t')-^- f [g{t'),c\t')] + dt> . (132) 

Jt Jt„ 1 Jt 

Notice that stochastic integrals use the full operator g, while the correction terms all use the 
restricted system operators g, c. 
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Bosonic integrand g: In this case the restricted operator g is identical with the full system 
operator g, so we can write everything in terms of g. 

(S) f g(f)dF(t') =00 f g {t')dF{t') + ^- f ' \g{t'),c{t')]dt' , (133) 

Jto Jto 1 Jto 

(S) [ dF(t')g(t') = (I) [ dF(t')g(t')- ^ {1 ~ N) f \g(t'),c(t')]dt' , (134) 

(S) f g{t')dF\t')={\) f g{t')dF\t')+ ^ {1 ~ N) f [g{t'\c\t')]dt' , (135) 

Jt Jto 1 Jto 

(S) f t dFUt')g(t') = a) [ dFUt')g(t')-y^- [ \g(t'),cUt')]dt> '. (136) 



to Jto Jto 



4.4. Stratonovich quantum stochastic differential equation 



We want to write the Stratonovich equivalent of the quantum stochastic differential equations 
( II 1411 15t . from which we see that the Ito integrals to be converted into Stratonovich integrals 
correspond to making the choices 

g(t)dF(t) -> -^[a{t),c\t)]I{t)dF{t), (137) 
dirt (t) $(*)-> SfdF\t)I(t)[a(t),c(t)]. (138) 

The combinations [a(t),c^(t)] i(t) and I(t) [a(t),c(t)] thus form the appropriate Fermionic 
full system operators. Using therefore dl29ll32| i. we find the total correction term is 

-^2{[c,[a,c^}] + -[c\[a,c}] + }, (139) 



and this is precisely the negative of the third line of il 15i . which leads to the Stratonovich 
form 

(S)da = -^[a,H sys ]dt+ ^ [a, c] - l[a, c f ]c} dt 

- V7 [a, c f ] IdF(t) + ^fdF\t) I [a, c]. (140) 

This corresponds exactly to the Langevin equation for the restricted system operator as given 
in (1721 . and thus justifies the form (II 1411 1 5i chosen for the Ito quantum stochastic differential 
equation. 

5. The master equation for the system density operator 

In the Heisenberg picture, the density operator is time independent, and takes the form 

P = Psys(to) <S> PF ■ (141) 

This corresponds to an assumption that the initial density operator in the Schrodinger picture 
can be factorized into a bath and a system term. 
Using this density operator, we would say that 

(a(t)) = Tr sys {ap sys (i)}, (142) 
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On the other hand, the quantum stochastic differential equation M 1 4I > can be used in 
conjunction with the fact that the means of the Ito integrals are zero to show that 

d_ 

dt 



-(a(t)) = (^-^{a,H sys \ + ^(l-N){2Jac-aJc-Jca} 



+ {2cac f -acct -cJa}j (144) 
from which we deduce, since this equation is true for any restricted system operator a 

= ^[/°sys, H sys ] + ^(1 - 7V){2cp sys c t - p sys c f c - c f cp sys } 

+ ^■{2c t /9sy S C - cS p sys - PsysCC 1 }. (145) 

6. Quantum stochastic differential equations in the interaction picture 

For a discussion of many aspects of input-output theory it is advantageous to formulate 
an appropriate quantum stochastic differential equation theory in the interaction picture, as 
explained in Chap. 1 1 of ]5) in the case of Bosonic noise. In such a situation, the time- 
dependent state vectors can be written in terms of the evolution operator U(t,to) as 

\1>,t)=U(t,to)\il>,to)- (146) 

In much the same way as for Bosonic noise, we can write a Stratonovich quantum stochastic 
differential equation for the evolution operator corresponding to the Hamiltonian OH in the 
form 

(S)dU(t,t ) = -^H sys U(t,t )dt+{^dFi(t)c + ^fdF(t)<*}U(t,t ). (147) 

Our aim now is to transform this to a corresponding Ito form. We assume there exists an Ito 
equation in the form 

a)dU(t,t ) = {a(t)dt + P(t)dF f (t) + ft(t)dF(t)}U(t,t ), (148) 

and from this derive the relationship between the Stratonovich and Ito integrals. 
A Stratonovich integral of the evolution operator is defined by 

(S) f dF(s) U(s) = lim £ m U(U + ,) + U(U) _ 

J n—*oo z — J Z 

We now express Ui+i in terms of Ui using the Ito quantum stochastic differential equation 
dl48> . and neglecting terms of order of magnitude At^ 2 and higher, we get 

(S) J dF(s) U(s) = Jim AF * j 1 + \^ AF i) ^(*0 ( 150 > 

= lim V { AFi - - N)Au] U{U) (151) 

= (I) / dF(s) U(s) - hi — N) f (3(s)U(s) ds. (152) 



2 

Similarly 

(S) J dF*(s)U(a)=Q) J dF^(s) U(s) — J ft (s)U(s) ds. (153) 
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Using these relations, we can then convert the Stratonovich quantum stochastic differential 
equation to the equivalent Ito form 

(X)dU(t,t ) = - Un sys + ^22^ + 7(1 2 ^ gtgj U(t,t )dt 

+ {^dF\t)c + ^jdF{t)$} U(t,t ), (154) 
The full density operator at time t can be written 

p(t) = U(t,t )p(t )U^t,t ) (155) 
and so obeys the quantum stochastic differential equation 

dp(t) = {-L[H sys ,p(t)} - l^[ P {t)^] + l{1 ~ N) [p(t),c^] + \dt 

- -fcdF^(t)p{t)dF(t)^ - jcUF(t)p(t)dF^{t)c 

+ yft[dF\t)c + dF(t)&,p(t)] (156) 

6.1. Alternative derivation of the master equation 

To derive the master equation from the evolution equation dl56> still requires us to use the 
restricted system operators as follows. Firstly, note that the for the product operators we have 

c^c = c'c, cc' = cc\ (157) 

and of course the restricted operator form of H sys is H sys itself. These operators are therefore 
proportional to the identity in the bath space, and we can write, for example 

Tr B {[p(t),cc*}+} = Tr B {[p(t),cc+] + } = [Tr B {p(t)}, a?]+ = [p sys (t), cc*]+ (158) 

However, for the terms on the second line of dl56> this is not immediately possible, since c 
and do act in the bath space, since they do not commute with the noises. We therefore have 
to convert to restricted system operators . In the interaction picture the operator / introduced 
in J64i is time independent, so we can write, for example 

Tt B {cdF*(t)p{t)dF{t)&} = Tr B {cl dF^ (t)p(t)dF(t) 7c f } 

= cTrs {I dF\t)p(t)dF{t) 1} c ] 

= cTr B {dF(t) I 2 dF^(t)p{t)} c ] 

= (l-N)cp sy J. (159) 

The trace over the last line of dl 561 gives zero, so we derive again the master equation in the 
form 11451 . 

6.2. Correlation functions 

From the density operator one can in principle compute all correlation functions for the 
restricted system operators. Using the relationships between the operator forms a, a and 
a, as applied in Sect l3.3.41 it is then possible compute the correlation functions for the full 
system operators. 

The importance of the full system operators comes from the need to compute the 
correlation functions of the output operators which arises because of the relationship ( I59> 
between outputs, inputs and the full system operators. 
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6.2.7. Vacuum input We can compute 

OoUO/outW) = {{flit') ^jc\t')){Ut) V7c(t))) • (160) 

The simplest case is if the input field corresponds to the vacuum, in which case we get 

(fl t {t')foAt))=l^{t')c{t)) (161) 
= 1 (c\t')c{t)) (162) 

If we want the number counting correlation function, then we find in much the same way as 
the Bosonic case that for the time ordered correlation function, in which t' > t 

(/oU*)/oUO/out(0/out(*')> = 7 2 (5 f (t)c\t')c(t')c(t)) (163) 

= 7 2 {<*(t)c\t')c(t')c(t)). (164) 

6.2.2. White noise input Let us consider again the correlation function (/o U t(0/out(£)) in 
the case that the input is non-vacuum, and for t < t'. In this case (I160> becomes 

(fL(t')font(t)) = ((4(0 - ^(OK/inW - y/im) d65) 

= -V7<5 + (0(/mW-V7c(t)) (166) 
= - V7 <^(*') W*) - n/7c(*) dt}) /<fc. (167) 
Here we will consider dF(t) to be an Ito increment, but there will be no difference between 
an Ito and a Stratonovich version of the increment in this formula except when the two times 
t and t' are equal. 

Thus it is necessary to compute {& (t 1 ) dF(t)) where t < t'. For a general full system 
operator a(t') we can write (where t\ = t + At, for brevity) 

(~a(t')AF(t)) = 

Tr sys {Tr (tlit ,) {Tr (Ml) {aU(t\ ti)U(ti,t)AF(t)p(t)U~ 1 (ti,t)U~ 1 (t' 7 ti)}}} 

(168) 

The evolution operator U(t' , t\) contains no dependence on the noise in the interval (t, t\), so 
we can than write 

(5(0 AF(t)) = 

^{aTr^^^jt/^^OTr^^^l^^^OAFW^f/- 1 ^^)}^- 1 ^,^)}} 

(169) 

We now write the infinitesimal form 

U(h , t) = U(t + At, t) = - | ^H sys + + 7(1 2 ^ ^gj A< 

+ {V7AFt(t)5+V7AF(t)5t} (170) 

Now computing the trace over the interval (t, t\) = (t,t + At), and keep only terms of order 
At, we get 

Tr (tiil) {C/(t 1 ,t)AF(t)p(t)[/- 1 (t 1 ,t)} =^yN[c,p(t)]At, (171) 
so that 

(d(t')AF(t)) ^^NTv sys {hTT {tun {U(t\t 1 )[d,p(t)}U-\t',t 1 )}}At (172) 
Taking account of the fact that for t' < t the average obviously vanishes, we get 

(a(t')dF(t)) = ^Nu(t' -t){[a(t'),c(t)]±}dt. (173) 
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Using this result we get 

(/oUO/outW) = 7(1 - N)(cHt')c(t)) - yN{c(tycHt')) + N5(t t'). (174) 

The result is valid for all t, t'\ the case for t > t' is given by the complex conjugate of that for 
t < t'. 



7. Examples 

7.1. The two-level ion 

We consider a two level system in which we have a lower energy level with an even number of 
electrons and an upper level with an odd number of electrons — hence the terminology "ion". 
The system interacts with an electron field, and thus may be described by the choices 

c = a~, (175) 
ct =er+, (176) 

~2 



H sys = — <t z , (177) 



so that 



c = Ia~ = a-, (178) 

c f = Ia+ = a+, (179) 

J=<Tz, (180) 

c=J<j-=-a-, (c)t = a+J = (181) 



Quantum Langevin equations for the full system operators Using the usual commutation 
relations in J53i . we get the explicit equations of motion 

(iw - |) - V7f?M d82) 

&~ = - (iw + |) cr- - V7/inW, (183) 

o z = - j(a z + 1) - 2^a+f in (t) - 2^ a- fl{t). (184) 

The equations Jl 8211 83^ are linear in a ± , f\ n , fl v and are thus exactly solvable. Although 
the equation for a z is not linear, it solution follows from the other equations by using identity 

[a+,a-] = a z . 

This solvability arises because a two level ion is in fact the same thing as another Fermion 
degree of freedom, that is, in this case the operators c, & are like Fermion creation and 
destruction operators since [c, c^] = 1 and they anticommute with the bath Fermion operators. 
This is exactly the same situation as for a harmonic oscillator interacting with a Bose input 
field, which is exactly solvable in the same way. 



Stratonovich quantum stochastic differential equations for the restricted system operators 
These take the form, from (11401 

(S)da+= (\io-^)a+dt + ^dF\t){Ia z ), (185) 

(S)da~= - {vjj+T^a- dt + ^ (Ia z ) dF(t), (186) 
(S)da z = -j(a z + l)dt + 2^o- + (Ia z )dF(t) + 2^dF i <(t)(I<T z )<j-. (187) 
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Note that from Sect l3.3.41 the quantity which multiplies all the noise terms, Icr z is in this case 
the conserved operator K. Thus the first two equations are again exactly solvable. 

It is surprising to see the apparent difference between the form of the set of equations 
( 1182141841 and that of the set H85H187i — however, it must be borne in mind that the 
transformation from the full system operators to the restricted system operators involves the 
time dependent operator I(t), and this accounts for the apparent difference. Moreover, once 
the conserved nature of K = Ia z is noted, one sees that the only difference between the 
two sets of equations is a sign change when the total number of Fermions is odd, which is of 
course of no consequence. 

Ito quantum stochastic differential equations for the restricted system operators These 
take the form, from J 1 1 51 . 

(I)da + = (iio-^jn+dt + ^dF\t){Ia z ), (188) 

(I) da- =-(hj+ |) a- dt + V7 {!*,) dF(t), (189) 
(I)da z = - 1 (a z + l-2N)dt + 2 y /j<T + {I<T z )dF{t) + 2^dF\t)(Ia z )<j-. (190) 



7.2. The harmonic oscillator 

The harmonic oscillator coupled to a fermion bath is a problem which cannot be solved 
exactly, in the same way as the two level atom coupled to a Bosonic bath is not exactly 
solvable. In this case we have harmonic oscillator creation and destruction operators a, a', 

c =a, (191) 
C* =a\ (192) 
H sys = huala, (193) 



so that 



c =Ia = a, (194) 

& = la* = 5 f , (195) 

J = (-l) a1a , (196) 

a = Ja, (a)t = a f J. (197) 



Quantum Langevin equations for the full system operators Using the usual commutation 
relations in ( I53l l, we get the explicit equation of motion 

5= -i^5-(2a t a+l){|a + V7/m(t)}-{^5 t + x /7/ i t n W}(2a 2 ). (198) 

Stratonovich quantum stochastic differential equations for the restricted system operators 
We get the form, from (11401 

(S)da= - (icj + adt- yfyIdF(t) (199) 

The difference between the equations for the full and the restricted operators is now quite 
dramatic. The simple form of dl99l is deceptive, since it involves the time dependent operator 
/, whose equation of motion is not solvable. 
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Ito quantum stochastic differential equations for the restricted system operators We get, 
from fTTBV 

(I) da = - (iw + 7( ' 1 ~ 2jV ' ) ) a dt - ^ I dF(t) (200) 



2 

(l)d(a t a) = - 1 {(l-2N)a t a + N}dt- y /ja^IdF(t)-^/jadF t (t)I (201) 
Using the last equation, we can find the equation of motion for the mean number n(t) to be 

n(t)= - ^{(1 - 2N)n(t) + N} (202) 
with the stationary solution 

N 

(203) 



1 - 2N 

For Fermionic noise at a predominant frequency u>, we know that 

N = , lur (204) 

and this yields the correct result for the harmonic oscillator 

1 



s e hcj/kT _ i ' 

8. Conclusions 



(205) 



This paper answers what is in some sense an academic exercise — how do we deal with the 
Fermion inputs and outputs that are so common in the real world, turning up in electronic and 
many other systems. The two previous treatments 0E1 gave partial answers, the first dealing 
only with counting statistics, the second dealing only with noninteracting particles in cavities. 
The treatment given here is compatible with both. 

As coherent Fermion physics becomes more important, for example in highly degenerate 
trapped cold Fermi vapors |6| this kind of formalism will undoubtedly become more relevant. 
I look forward to that time. 
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